Numerical Study on Spontaneous Symmetry Breaking in a XY Quantum 
Antiferromagnet on a Finite Triangular Lattice 
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Motivated by recent experiments that require more complicated macroscopic wave functions in 
the condensed matters, we make numerical study on a XY quantum antiferromagnet on a finite tri- 
angular lattice using the variational Monte Carlo method and the stochastic state selection method. 
One of our purpose is a numerical confirmation on dominance of a Nambu-Goldstone boson in low 
energy excitation. For another purpose, we calculate energy, an expectation value of a symmetry 
breaking operator and structure functions of spin by fixing a quantum number of the symmetry. 
These calculations are made for states that become degenerate in an infinitely large lattice. 

By numerical calculations we confirm existence of a Nambu-Goldstone boson, and find dependence 
of a square of the quantum number for the above quantities. Using these results we can discuss on 
complicated macroscopic wave functions in quantum spin systems. 



I. INTRODUCTION 

It is well known that the ground states of many quan- 
tum antifcrromagnets on two dimensional lattices exhibit 
semi-classical Neel orderl, which has been supported 
strongly by the spin wave theory(SWT)^ as well as nu- 
merical works^. This order implies that the ground 
state is a coherent state that consists of highly degen- 
erated states. Also this ground state is characterized by 
order parameters, which correspond with a magnitude 
and a phase of a macroscopic wave function. 

However, experimental works in other condensed mat- 
ters give us more complicated phenomena. For exam- 
ples experiments with alkali atoms have realized that two 
or more Bose-Einstein condensates with different phases 
merged and produced an interference pattern in their 
densities^. This interference pattern has forced us to 
examine again theoretical descriptions based on the co- 
herent stated""—. A work on superconductors^ can be ref- 
ereed as another experiment on the interference between 
macroscopic wave functions. By these experiments we 
have to recognize that the ground state or the macro- 
scopic wave function in the condensed matters, where 
spontaneous symmetry breaking(SSB) occurs, can be not 
described only by a few parameters. 

Although above experiments have not been done yet 
for quantum spin systems, at least many-body systems, 
theoretical investigations are needed for the ground states 
in these systems if we apply more complicated external 
interactions. For these investigations we focus our study 
on an effect due to finiteness of a system and a Nambu- 
Goldstone(NG) boson. In a finite lattice we does not have 
degenerate states, but have a state with the lowest energy 
for a fixed quantum number of a continuous symmetry. 
We would like to make numerical study for this state 
with the fixed quantum number. By this study we can 
make a theoretical discussion on the ground state with 
the complicated interaction. Also cluster experiments of 
quite small sizesi^ are another motivation of this study 
on a finite system. 



Although we have calculated energy of the NG boson 
by the SWT^, we have to assume that the SSB occurs 
for an application of the SWT. While in numerical works 
such as a Monte Carlo method, we do not need this as- 
sumption. But a confirmation of a NG boson in these 
works is not an easy task because a NG boson exists 
only for a small wave vector whose calculation requires 
a quite large lattice. Our numerical study will be made 
on a 108 site lattice, where the smallest magnitude of 
wave vectors is 27r/9, so that we can present quantitative 
discussions on the SSB of quantum spin systems. 

We calculate energy, an expectation value of a sym- 
metry breaking operator and structure functions of spin, 
which give us knowledge on the NG boson. In order to 
make clear our calculations, we introduce notations. If a 
Hamiltonian H in the system has a continuous symmetry, 
we define a charge operator Q whose quantum number 
can be fixed for the eigen state of the H . We define a 
state I n) for a system with a site number N . 



H I n) 



n)E{N, n) 



Q \ n) —\ n)n . 

Here E{N, n) is the lowest eigenvalue for a fixed charge 
number n. These | ri)'s become degenerate when a lattice 
size is infinitely large. Also in a finite lattice for many 
antiferromagnet systems, one could not see degenerated 
states for the lowest energy so that we have only one state 
of the lowest energy. When we add an external operator 
B to the H, by a modified Hamiltonian Hb = H + B, 
we obtain an eigen state | E{N, B)) of the lowest energy 
E{N,B) of Hb- Here | E{B)) = E„ | n)c„(B). If we 
control a form of B and its magnitude, we have the state 
of I E{N,B)) with various coefficients Cn{B). Cn{B) is 
determined by E{N, n) and {n \ B \ n') . 

In usual experiments it is difficult to control B so 
that we assume that B = h'^^[aj^{i) + &-{%)), where 
{n + 1 \ a+{i) \ n) — V ^ and (T-{i) — so that 

we have the coherent ground state where the coefficient 
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Cn{B) is Cexp{—in9). However, recent experiments^ 
have told us that one can control B . Especially in a 
superconducting single-electron transistor— the ground 
state has only a few | ri)'s. These experiments ask us 
a following question; What kind of B induces a compli- 
cated ground state which shows phenomena differed from 
that found in the coherent ground state? A definite an- 
swer to this question is a final goal of our study. The first 
step to this goal is to examine extensively n-dependence 
of expectation values of operators for | n) . 

In the finite lattice system with Hb, an expectation 
value of an operator O is given by 

{E{N,B) I I EiN.B)) = ^(n | 6 \ n')c„(B)c„- (B) . 

n,n' 

Here we notice that {E{N, B) \ 6 \ E{N, B)) consists 
of a n-independent term (O) and n-dependent terms of 
A{0){n,n') of (n I O I n'). For the coherent states it is 
assumed implicitly that we observe only (O). For small 
N systems or complicated macroscopic wave functions, 
it is possible to observe A{0){n, n'). Therefore we would 
like to study on (n | O | n'). These calculations have 
not been made yet, at least to my knowledge, except the 
energy. For the energy, previous works have confirmed 
the n^/N dependence numerically^ i ^ ^ i ^ ^ . 

Our study is made on the XY quantum antiferromag- 
net on a triangular lattice because of foUowings reasons. 
The first reason is that this system has the U(l) sym- 
metry, the simplest one among continuous symmetry 
groups, which are indispensable for existence of the NG 
boson. Due to the same reason there has been many 
works on the XY quantum antiferromagnet, specially 
the square lattice modeU^ii^. They support that the 
XY model gives us essential properties of the Heisenberg 
model. 

Another reason is that previous works have found the 
excellent trial state for the XY antiferromagnet on the 
triangle lattic&l^'-^''. We can expect that the variational 
Monte Carlo(VMC) method is powerful in this system. 
By the VMC method it becomes possible to study the NG 
boson in a large lattice, which is needed for study of exci- 
tation with small wave vectors. Since the trial state is es- 
sential, we must examine reliability of this state. For this 
examination we use the stochastic state selection(SSS) 
mcthodi^i"— . Note that it is difficult to calculate only by 
the VMC method. 

A plan of this paper is as follows. In section 2 we de- 
scribe the model and calculation methods. In section 3 
we present numerical results. Lattice sizes are TV = 48 
and N ^ 108. Calculations in TV = 36 and = 324 
are made for energy. After fixing parameters of the trail 
state in the first subsection, we show the lowest energy 
for the state with a fixed value of the charge Q in sub- 
section 3.2. Also we employ the SSS method to estimate 
quality of the trial state used in study by the VMC. In 
subsection 3.3 we present results for expectation values 
of a symmetry breaking operator. In a next subsection 



we will show results of the structure functions, which are 
determined by the property of the NG boson'^'^. In sub- 
section 3.5 we will give a direct calculation of energy of 
excitation at a small wave vector. By this calculation we 
obtain a velocity in this system^i^— . A final section is 
devoted to a summary and discussions. In appendix A 
we make a brief description on the VMC method and the 
SSS method. Here we describe a way to connect the SSS 
method with the VMC method. Another appendix gives 
us a discussion on n-dependence of the expectation value 
of the spin operator. 

II. MODEL AND CALCULATION METHODS 

A system which we study is the quantum XY antifer- 
romagnet of spin one-half on the triangular lattice. Its 
Hamiltonian is given by 

H = J2iSA^)SAJ) + Sy{^)SyU)) ■ (1) 

Here Sxii), Sy(i), Sz{i) are x-component, y-component 
and z-component of one-half spin operators on the i-th 
site and the sum runs over all bonds of the A^-site lattice. 
This Hamiltonian has the U(l) symmetry, whose charge 
operator is defined by 

Q= J2 ^^w- (2) 

i=l,...,N 

Since we would like to make calculations in a large 
lattice, we use the VMC method. In the XY spin system 
on the triangular lattice, previous works give us the good 
trial state, | ^'). 

I - E I {^a>c({*^}) ' (3) 

where | {si}) is a basis state. Here {si} — (si, S2, • ■ • , sn) 
, where Si = —1/2, or 1/2 . The trial state^^ is given by 

I *) = E Texp{glJ29^3S^SJ) I {sj) , (4) 

(si) id 

T3 = exp{il3 E lijkSiSySk) , 

i.j.k 

gij = E ea;p(ifc • {xi - Xj))v{k) , 

k 

v{k) = 1 - ^ , 
Vl-7(fc) 
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= {cos{k^) + 2cos{k^/2)cos{V3ky/2)}/3 , 

where we should note that each site is categorized 
into three sublattices, A-sublattice, B-sublattice and C- 
sublattice. Also jijkSiSjSk is a three body interaction, 
which is given in Ref.^^. Note that | is a sum of | n)^. 

I *) J2 I ")*^« ■ 

n=-N,...,N 

Here Q | n)^ = n | . This trial state has two param- 
eters, gi and /3. We search the minimum energy state of 
n = by changing these parameters. After fixing val- 
ues of gi and /3, we calculate an expectation value of the 
Hamiltonian squared in order to estimate quality of the 
trial state. For these calculations we use the SSS method. 
Then we will confirm that a difference between the expec- 
tation value of the Hamiltonian squared and the square 
of the expectation value of the Hamiltonian is small. 

In the infinitely large lattice, properties of the NG bo- 
son come from a following equation 

{G\[Qj{i)]\G} = {G\S^{i)\G} , 

using the ground state \ G). If the right-hand expectation 
value is not zero, the SSB occurs. However for a finite 
size system, we have only one state \ n = 0) for the lowest 
energy and this state is the eigen state of Q. Therefore 
this expectation value vanish for a finite size lattice. In- 
stead of the equation, we use a following equation for a 
finite size lattice. 

(m| [Qj{i)] I n) = (m I 5^(0 I n) . 

When the right-hand side of the above is not be zero, 
we see a signal of the SSB in a finite size lattice. In our 
study we adopt S+{i) — Sx{i)+iSy{i) as The above 
equation becomes 

{n + l\[Q,S+{{)]\n) = {n + l\S+{t)\n} . (5) 

This trial state | ^P) is constructed using an assump- 
tion that each site is categorized into three sublattices, 
A-sublattice, B-sublattice and C-sublattice, where direc- 
tions of a unit vector are 0, 27r/3 and 47r/3 in a; — y plane 
for these sublattices. Therefore the expectation values of 
S+{i) for the i-th site on a sublattice differ from one on 
the another. For getting the same value for every site, 
we make a rotation around the 2-direction of spins. For 
the i-th site on the A-sublattice we make no rotation. 

= SA^) , S^{^) = Sy{^) ■ 

For the i-th site on the B-sublattice we make a rotation 
with an angle 27r/3. 



For the i-th site on the C-sublattice we make a rotation 
with an angle 47r/3. 

S^i^) = ~lsA^)+^Sy{^) , = -ylsA^)-^Syi^) . 

In this representation we have the same expectation value 
of spin operators, S!^{i) = S^{i) -\- iSy{i) , being inde- 
pendent of the site. 

{n + l\S'^{i)\n) =v{n) . (6) 

Note that this expectation value depends on n. 

Also using these the rotated spin operators S^{i), 
Sy{i), S^{i), we calculate structure functions. For this 
purpose we introduce spin operators that depend on wave 
vectors. 

S^ik) = E e^Pi^kn)S^{z) . (7) 

i 

Here fc is a wave vector, fi — ni(i)ei -I- n2{i)e2 is a site 
vector, which is defined by integer numbers ni(i) and 
n2{i), and r = x,y,z. Also ei = (1,0), 62 = (l/2,-\/3/2) 
are unit vectors on the triangular lattice. In future de- 
scriptions of S^{k), we omit the superscript R in order 
to avoid complexity. 

In order to make clear a relation between the NG boson 
and the structure function, we define the NG boson. A 
boson operator is given by the annihilation a(fc) and 
creation operators a'{k). 

<l(i) = -^= — j^^^[exp{ikri)a{k)-\-exp{—ikfi)a^{k)]. 
^2E{k) 

(8) 

where E{k) is a energy, and in the infinitely large size 
lattice, E{k) — c \ k \ for a small | A; |. As previous 
works show'^^, the NG boson appears in the charge cur- 
rent Sz{i)- 

Szii) = fdt^{i) -\- other terms . 

Here dt is a time derivative and we can neglect 
other terms for a small wave vector. Also / is a con- 
stant. From these discussions we have an equation in the 
infinitely large lattice. 

(G I SMS.{k) I G) = ^^^^ , (9) 

for a small | A; | . We will examine this equation carefully 
in the finite size lattice. 

Also the field theoretical argument shows that the NG- 
boson appears in the spin operator S'y(i), 

Sy{i) — Z^{i) -I- other terms . 
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N 


36 


48 


108 


324 


91 


0.07 


0.05 


0.023 


0.0075 


E{N,0)/N 


-0.40794 
±0.1 X 10""* 


-0.406718 
±0.8 X lO"'^ 


-0.405383 
±0.4 X lO"'^ 


-0.404905 
±0.15 X 10-"*" 



B. Energy 



TABLE I. Values of parameter gi of the trial state and the 
minimum energy for the lattice size A''. 



Therefore we have 



(G I Sy{k)Syik) I G) 



2c I k 



(10) 



for a small | k \ and the infinitely large size lattice. We 
will study the above equation in subsection 3.4. 

III. RESULTS 

A. The trial state 

As described in section 2, the trail state is determined 
completely by parameters gi and /3. By changing values 
of gi and /?, we find the minimum value of the expectation 
of the Hamiltonian. By the minimum value we obtain the 
best values of these parameters. As said previously, the 
trail state is given without fixing values of Q. That is 



N/2 

E 

= -N/2 



n)'q,Cn 



In a following discussion we omit a subscription "^f" of 
a state | n) in order to avoid complexity. Here the state 
I n) is a complex state. 

I n) =1 n , R)+ \n ,1 )i . 

For finding the minimum value of the expectation value 
of the Hamiltonian, we use only a real state \ , R) 
for avoiding statistical fluctuations due to the redundant 
complexity. After finding the minimum value, we have 
confirmed a orthogonality of the real state and the imag- 
inary state, and the same value of the expectation value. 

(0,/|0,i?)==0 , 

{0,R\H\0,I)^{0,I\H\0,I) . 

The expectation value of the Hamiltonian is denoted by 

E{N, n) = {n ,R \ H \ n ,R) . 

As results we obtain /3 = 0.09, which is the same value for 
all lattices with N = 36, 48, 108, 324. Obtained values for 
gi are given in Table 1. If (71 is a function of the lattice 
size N, we guess a function form, 

gi{N)r^0.25/N . 



First we show results on the minimum expectation 
'alue of H for charge n = in lattices of various size 
V , which are given in Table ID We would like to examine 
reliability of these results. For A'' = 36 we can compare it 
with the exact value, Eq/N = —0.41095 that is obtained 
by the exact diagonarizationiS. A difference between the 
exact value and the VMC value is about 0.7% 

Next on examinations for N — 48, 108, we use the SSS 
method to calculate expectation values of the Hamilto- 
nian square H^. By this calculation, we obtain a follow- 
ing ratio. 

SiE) = i{0\H'\0)^i{0\H\0)f)/i(0\H\0)r . 

This ratio gives us rough estimations on differences be- 
tween results and the exact values. 

(5(£;) = 1.1 X 10"^±0.2 X 10"^ /or TV = 36, 



S{E) = 4.1 X lO"-* ± 1.6 X 10"'' 



for iV = 48 



= 4.2 X 10""* ± 0.4 X lO"'' f or N = 108. 

For = 36 the ratio 0.11% is somewhat smaller than 
the difference 0.7%. Quantitative estimations on the dif- 
ference from the exact value are difficult, but quite small 
values in these results justify our study on the SSB by 
the VMC. 

From obtained energy we can have size dependence of 
the energy, which is 

E{N,0)/N ^ ea+a/N^^^ , 

where eo = -0.4048 ± 0.0001 and a = -0.65 ± 0.02. 

Next we will show results on expectation values of H 
for non-zero values n. Previous studies based on the 
SWT and numerical approaches show that dependence 
of the energy on n is 



E{N, n)/N = E{N,Q)/N + br?/N'^ 



(11) 



For A^ = 48 and A^ = 108 we plot E{N,n)/N as in 
FigUl These results strongly support the above depen- 
dence on n. Also we should note that the above depen- 
dence pT|) is acceptable for quite large n. By the least 
square fitting, they are b = 1.715 ± 0.004 for A^ = 48 and 
b = 1.726 ± 0.002 for A^ = 108. 



C. Expectation values of S+{i) 

As described in section 2, an expectation value v{n) of 
a operator S+{i), which show the breaking of the U{1) 
summery, is given by Eq. ([S]) . This equation assumes that 
these are independent of sites and real. In order to verify 
these assumption, we calculate a standard deviation on 
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sites and imaginary parts of v{n). The site averages of 
real and imaginary parts of the expectation value for n — 
are 

(1 I S+i^) I 0) = 

0.4815 ± 0.0002 + i(-1.27± 1.40) X 10"'' /or = 48 , 
(1 I S+i^) I 0) = 

0.46425 ± 0.000008 + i(-l.ll ± 1.52) x 10"^ for N = 108 

From estimations of the statics error, we can say 
that the standard deviations on sites are 1.4 x 10^^ = 
2 X IQ-'^Vn for TV = 48, and 8.3 x 10"'' = 8 x lO^^ViV 
for = 108. These values show that the above assump- 
tions are justified numerically. Dependence of expecta- 
tion value v{n) on n is shown in Fig. [2] Here the hori- 
zontal axis is denoted by n(n + 1), and a vertical axis is 
done by {v{n)}^. We find that this dependence is well 
described by a linear function of n(ri -|- 1). 



[w(„)]2 ^ [v{0)]^+d 



n{n + 1) 

]V2 



Values of d are -0.82 ± 0.02 for N = AS and -0.75 ± 0.06 
for N = 108. As discussed in Appendix B, we have d = 
— 1 if we assume that we neglect contributions of the 
excited states when we calculate expectation values of 
[S'+(i), S-{j)] — 2iSijSz{j)- This small discrepancy on d 
shows that we should not neglect these contributions. 



D. Structure functions 

In this subsection we show results ori^ the structure 
function on a product of Sz{—k) and Sz(k), which is de- 
fined in a finite size system by 

Fzz{k,n)^{n\Sz{~k)Sz{k)\n) . (12) 

For I '$), we have 

(* I Sz{-k)Sz{k) I *) = ^(n I Sz{-k)Sz{k) I n)CZCn 

n 

= ^F,,(fc,n) I C„ |2 . 

n 

If the NG boson dominates in the structure function for 
small wave vectors, 

Fzz{k,n) — — c I k I +n — dependent terms . 

To check this form, we plot results on Fzz{k, n) as a func- 
tion of I fc I in Figures m IH For | k |< 2, we can see a 
linearity of | fc | for TV = 48 and iV = 108. This linear- 
ity strongly supports that the NG boson dominates in 
the low energy spectrum. Also we find that Fzz{k, n) de- 
pends on n. In order to see this n-dependence in detail. 



we plot Fzz{k,n) as a function of n^. From these plots 
we find that Fzz{k,n) are linear functions of n?, which 
are 

Fzz{k,n)^Fzz{k,0)+gzz{k)n^ ■ 

Next we plot gzz{k) as a function of | fc |. By results we 
guess that this is a linear function of | fc |. We make the 
least square fitting on this function using data of | fc |< 2. 

Fzz{k,n) = i^c\k\+{azz + l3zz\k\)nyN^ . (13) 
From these figures we obtain for = 48 and N — 108, 
4rc = 0.130 ± 0.004 , a^^ = 0.26 ±0.07 , 



, = -0.60 ± 0.06 for TV = 48 , 

c = 0.126 ±0.003 , = 0.45 ±0.07 , 
^ = -0.67 ±0.08 for TV = 108 . 
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A squared error of the expression p3)) is calculated by 
an average of square of differences between data and the 
expression. The error par one data is 0.0036 for TV = 48, 
and 0.0025 for TV = 108. These small values imply that 
the function ([T3)) can describe quite well the data. 

Next we show results on the structure function on a 
product of Sy{—k) and Sy{k), which is defined in a finite 
size system by 

Fyyik,n) = -{n + l \S+{-k)S+ik) |n-l)/2 
- (n - 1 I S-{-k)S-{k) I n + l)/2 
+ (n I S-i-k)S+{k) I n)/2 

+ {n I 5+(-fc)5_(fc) I n)/2 . (14) 
For I ^P), we have 

(* I Sy{~k)Sy{k) I *) = 

Y.{-{n + 1 I ^+(-fcl5+(fc) I n - l)/2Q+iC„_i 

n 

- (n - 1 I ^_(-fc)5_(fc) I n + l)/2C:_iC„+i 
+ {n I 5_(-fc)^+(fc) I n)/2C*_,Cn+i 

+ {n I 5+(-fc)5_(fc) I n)/2QC„} 

- ^ Pyyik, n)C*Cn 

n 

- + 1 I S+{-k)S+{k) I n - l)/2(C*+iC„_i - Q^) 

n 

+ (n - 1 I ^_(-fc)5_(fc) I n + l)/2(Q+iC„_i - QC„)I15) 

This equation shows that if we know coefficients C„, we 
calculate the expectation of (^P | Sy{~k)Sy{k) \ 5'). 

If contributions form the NG boson dominate in the 
structure function for small wave vectors, 



Fyy{k,n) 



2c I fc 



n — dependent terms 
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To check this form, we plot results on Fyy{k, n) as a func- 
tion of 1/ I /c I in Figures m m For | fc |< 2, we can see 
a linearity of 1/ | fc | for iV = 48 and iV = 108. This lin- 
earity strongly supports that the NG boson dominates in 
the low energy spectrum. In order to see this n- depen- 
dence, we make the same discussion as that for Fzz{k, n). 
First wc plot Fyy(k,n) as a function of in?. Then from 
these plots we find that Fyy{k,n) is a linear function of 
n^, which are 

Fyy{k,n) ^ Fyy{k,Q) + gyy{k)n^ . 

Next we plot gyy{k) as a function of 1/ | fc |. Results 
suggest that this function is a linear function of 1/ | fc | . 
If we make the least square fitting on this function, we 
obtain 

Fyyi^^^)^T^ + iavv+fiyvl\^\)^VN^ ■ (16) 
2c I fc I 

From figures we obtain for = 48 and = 108, 

Z^/2c = 0.458 ± 0.004, Uyy = 0.90 ± 0.04, 
Pyy = -1.66 ± 0.06 for = 48 . 

Z^/2c = 0.439 ± 0.005, ayy = 0.89 ± 0.12, 
Pyy = -1.65 ± 0.15 for TV = 108 . 

A squared error of the expression (jl6p is calculated by 
an average of square of differences between data and the 
expression. The error par one data is 0.0027 for A^ — 48, 
and 0.0037 for A^ = 108. By these small values it is 
justified that the function p6|) describes the data as well 
as that for F^^ik, n). 

We make comments on a structure function of a prod- 
uct of Sy{—k) and Sz{k). In descriptions by the NG 
boson we have 

I Sy{-k)Sz{k) I *) =iZ//2 . 
Also from a commutation of 

[^^^(fci), Sz{k2)] = -j=Sx{ki + fc2), 
V A* 

we have Z f = v. In order to examine this equation, we 
calculate the structure function in the finite size lattice, 
which is defined by 

Fyz{k,n) = -J=(n + 1 I S+{k)Sz{~k) \ n) 

--^(n|5-(fc)5.(-fc)|n + l> . 

We have confirmed numerically that Fyz{k,nys are in- 
dependent of wave vectors fc and their values agree with 
v{n) within error, although results are omitted. 



E. Energy for excitations with small wave vectors 

In calculations of the structure functions, we do not 
see the energy uj{k) of the excitation, so that we could 
not determine a value of the velocity c, which is defined 
by c = a;(fc)/ | fc |. In order to obtain the velocity, we 
calculate the energy by a following method. If | fc, n) is a 
state of one NG boson with a wave vector fc and a charge 
n, we have H \ k,n) =\ k,n){u}{k) + E{N,n)}. Note 
that E{N,n) is the lowest energy for the charge n. It 
is difficult to get an eigenvalue so that we are contented 
with an expectation value of the H, which is given by 

cj^(fc) ^ {k,n \ H \k,n) - E{N, n) . 

The results on the structure function imply that the spin 
operator S'z(fc) can create a state of one NG-boson with 
a wave vector fc from the lowest energy state | n) . There- 
fore we can approximate the state of one NG boson by 
the lowest energy state that is operated by Sz{k). 

Sz{k) I n) -I k,n)A . 

Here A is a constant. If we use this approximation, a 
calculated energy iOA{k) is 

ujA{k)= {k,n\ H \k,n) - E{n) , 

= A^{k,n\H \ k,n)/{A^{k,n\k,n)) - E{n) , 
^ {n I Sz{~k)HSz{k) I n)/{{n \ Sz{^k)Sz{k) \ n)) 
- E{n) ■ (17) 

Results of calculations of w^(fc) are collected in Table 2. 
For n — the velocity is a constant for the small wave 
vector within some error. This value is a little larger than 
0.75 that is calculated by the linear spin wave theory. 
On the n-dependence of the velocity we could not make 
a definite conclusion due to the large error. Here we 
would like stress that these results are quite non-trivial 
because calculations on the Hamiltonian are completely 
independent from those on the structure functions. 

IV. SUMMARY AND DISCUSSIONS 

Recent experiments on condensed matters as alkali 
atoms and superconductors require more complicated 
macroscopic wave functions, which differ essentially from 
the simple coherent state. In order to understand phe- 
nomena by these wave functions in quantum spin sys- 
tems, we have made theoretical study on states of the 
SSB in these systems. In this work we used the VMC 
method for numerical study. In this method we have 
to assume that the trial state is close to an exact eigen 
state of the Hamiltonian. Based on this assumption a 
study was made for the XY quantum antiferromagnet 
on the triangular lattice, where the good trail state has 
been known. Also in order to justify this assumption, we 
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N 


n 


(kx, ky) 




c 


48 





(7r/3, 0) 


0.83 ± 0.04 


0.79 ± 0.04 


48 







1.46 ± 0.03 


0.80 ± 0.03 


48 


6 


(tt/S, 0) 


0.85 ± 0.02 


0.81 ± 0.02 


48 


6 




1 or* 1 r\ no 

1.3d ± 0.03 


0.75 ± 0.02 


48 


10 


(7r/3, 0) 


0.69 ± 0.04 


0.69 ± 0.04 


48 


10 


fin \ 


1.33 ± 0.03 


0.73 ± 0.03 


108 







0.56 ± 0.03 


0.80 ± 0.04 


108 







1.01 ±0.02 


0.83 ± 0.02 


108 







1.16 ±0.02 


0.83 ±0.02 


108 







1.40 ±0.02 


0.76 ±0.02 


108 


6 




0.58 ± 0.04 


0.83 ± 0.06 


108 


6 




1.02 ±0.03 


0.84 ± 0.03 


108 


6 


(4f-o) 


1.11 ±0.03 


0.80 ± 0.03 


108 


6 




1.42 ±0.03 


0.77 ±0.02 


108 


10 


(f ,0) 


0.53 ± 0.03 


0.76 ± 0.05 


108 


10 


*-3' 3V3' 


0.98 ± 0.03 


0.81 ± 0.03 


108 


10 




1.08 ± 0.03 


0.77 ± 0.03 


108 


10 


/ StT TT \ 


1.44 ±0.03 


0.78 ± 0.03 


108 


16 




0.52 ±0.2 


0.74 ± 0.04 


108 


16 


'>3'~3\/3'' 


0.94 ± 0.03 


0.78 ± 0.03 



TABLE II. energy of an excited state with a wave vector 



employed the SSS method to calculate the square of the 
Hamiltonian. By results on the square we confirmed that 
the trial state was of a high quality. 

Our numerical examination has been made on states 
that become degenerate in a infinitely large lattice, after 
the confirmation on dominance of a NG boson in low 
energy excitation. In the examination we calculated the 
energy, the expectation value of the spin operator and the 
structure functions of spin by fixing a quantum number, 
which is a value of the z component of all spins. 

Our results on numerical calculation in lattice sizes of 
48 and 108 showed that the energy of a small wave vector 
was linear to a magnitiidc of the wave vector. In addition 
it was shown that the expectation values of the operators 
studied here were linear functions of a square of the quan- 
tum number, . Though in our study this conclusion is 
limited to the XY model on the triangular lattice, we 
suppose that the same conclusion can be obtained for 
other quantum spin models because n-dependent terms 
are restricted partially by algebraic structures as seen in 
appendix B. 

Final comments are made on experimental observa- 
tions of the n-dependent terms. Our results show that 

for a ground state with n <^ N and a large N ^ one can 
neglect contributions of n-dependent terms. However if a 
external interaction is the charge operator B = hQ, the 
ground state is | n*), where n* = h/2d N. For this state 
there exists a possibility of observing the n-dependent 



terms. More detailed discussion will be made in future 
works. 



Appendix A: Appendix A 

In this appendix we discuss on the VMC method and 
the SSS method. A trial state | \E') is given by the coef- 
ficients c({sj}) on the basis state, as described in section 
2. 

In the VMC method, a probability variable X is a basis 
state of I {si}). A probability Pvmc{X = {sj}) for this 
variable is defined by 

PvMciX = {s,}) =1 c{{s,}) I /[E I c{{sj}) \] . 

In a Markov chain Monte Carlo method, what we need 
is only a ratio of probabilities, | c({s-}) | / | c{{si}) \ 
so that we do not need the normalization factor, which 
corresponds to the denominator. In this method a ex- 
pectation value of a operator O in the trial state | \E') is 
given by 

I 6 I I *) 

E{.s.} I c{{Si}) P 

= [ T.T.^^(i''i})0{s'Ms,}c{{sj})/ I c{{sj}) I 

■ PvMciX = {S,}) ] 

/[El ^({^^}) I Py^cix = {sj) ] , 

{si} 

where 

= iWi] I O I {sj}) . 

Here note that we do not need (\I> | \I>) = 1 . Actually 

the trial state we used in this study is not normalized. 

While in the SSS method we give a probability vari- 
able X^g.y to each basis state | {st}), whose value is 
or a^si}- Here l/a^st} = n7,in(l, | c{{si}) | C) < 1. 
A probability Psss{X[s,} = a{s.}) = > while 

Psss{X{s.} = 0) = 1 - An average of X^^.y 

is 1. We multiply this probability variable X^si} to a 
coefficient c({sj}) and we have X^g.-^c{{si}), that results 
in or phase{c{{si}))/C. Here phase{x) = x/ \ x \. 
A state with a nonzero coefficient is called as a sampled 
state after the multiplication. A value of C determines an 
average number of sampled states. That is, as C becomes 
large, this number increases and the statistical error de- 
creases. 

We will describe a calculation method for an expec- 
tation value of the Hamiltonian squared, {'^ \ \ '^). 
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First by the VMC method we make sampling to coUect 
Nsmpi basis states. Next we operate the Hamihonian 
H to these sampled states. After the operation, a num- 
ber Nh of the basis states with nonzero coefheients is 
a few hundred times of Nsmpi- If we operate H to Nh 
basis states, we need very huge number of basis states 
for calculations of its inner product with the trial state 
I . This huge number makes a calculation difficult due 
to limited resources of CPU time and computer mem- 
ory. In order to overcome this difficulty, we employ the 
SSS method to Nh basis states to reduce a number of 
basis states. After a reduction of the basis states, we 
operate H and calculate its inner product with the trial 
state I 4'). By repeating these calculations, we make a 
statistical average for (^J^ | \ 'i'). 

Appendix B: Appendix B 

In this appendix we will present a discussion on n de- 
pendence of v{n) = {n + 1 \ S+{i) \ n) . We start from a 
commutation relation on the spin operator. 

By the state | n) we calculate a expectation value of 
operators in both sides. 

(n|[^+(z),5_(j)]|n)=%(n|5.(z)|n) , 



(n I S+{i)S^{j) I n) = (n I S+{i) \ n-l)(n-l | S^{j) \ n) 



+ |n-l,e)(n-l,e|5_(j) |n) . 

e 

Here | n — 1, e) denotes an excited state of the charge of 
n—1. We assume that the second term can be neglected 
as described in Refi^. Noting that (n — 1 | S-{j) \ n) = 
{n I S+{j) I n — 1)''' = v{n — 1)\ we obtain a following 
relation by making a sum over i, j 



{v{n - 1)}^ - {v{n)y ^ 2n/N^ . 
Therefore we have 

{vin)f ^{v{0)f ~n{n + l)/N'^ . (Bl) 

Here note that v{n) is real as explained in section 3.3. 

The above discussion explains well the linear depen- 
dence of {v{n)}'^ on n{n + 1) which is found in Fig. [2l 
although decreasing rates differ from the above estima- 
tions even for N = 108. It may imply that it is too rough 
to neglect contributions from excited states in calculating 
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FIG. 1. The lowest energy for each n. Errors are small so 
that we do not plot error bars. The lines are given by applying 
the least square fit to data. 
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FIG. 2. Squares of the expectation values for each n. The 
bold lines are given by applying the least square fit to data. 
The slim lines are given by equation (|B1I) . 
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FIG. 3. Structure Function ^^^(fc) for iV = 48. The line is 
given by applying the least square fit to data of n = 0, for 
which the magnitude of wave vectors is less than 2. 
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FIG. 4. Structure Function F^^fc) for N = 108. The line 

is given by applying the least square fit to data of n = 0, for 
which the magnitude of wave vectors is less than 2. 




FIG. 5. Structure Function Fyy{k) for = 48. The hue is 
given by applying the least square fit to data of n = 0, for 
which the magnitude of wave vectors is less than 2. 
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FIG. 6. Structure Function Fyy{k) for = 48. The line is 

given by applying the least square fit to data of n = 0, for 
which the magnitude of wave vectors is less than 2. 



